Hole core in superconductors and the origin of the Spin Meissner effect 
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It is proposed that superconductors possess a hidden 'hole core' buried deep in the Fermi sea. 
The proposed hole core is a small region of the Brillouin zone (usually at the center of the zone) 
where the lowest energy states in the normal state reside. We propose that in the superconducting 
state these energy states become singly occupied with electrons of a definite spin helicity. In other 
words, that holes of a definite spin helicity condense from the top to the bottom of the band in 
the transition to superconductivity, and electrons of that spin helicity 'float' on top of the hole 
core, thus becoming highly mobile. The hole core has radius go = 1/2A_l, with Xl the London 
penetration depth, and the electrons expelled from the hole core give an excess negative charge 
density within a London penetration depth of the real space surface of the superconductor. The hole 
core explains the development of a spin current in the transition to superconductivity (Spin Meissner 
effect) and the associated negative charge expulsion from the interior of metals in the transition to 
superconductivity, effects we have proposed in earlier work to exist in all superconductors and to 
be at the root of the Meissner effect. 

PACS numbers: 



I. INTRODUCTION 



The Larmor diamagnetic susceptibility for carriers of 
density n in orbits of radius a is[l| 
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It is remarkable that the normal state Landau diamag- 
netism of metals is described by this expression for nor- 
mal state orbits of microscopic radius kp 1 : 
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with kp the Fermi wavevector, g{tF) — in/2tp the free 
electron density of states, and [lb — eh/2m e c the Bohr 
magneton. It is equally remarkable that the perfect dia- 
magnetism of the superconducting state is described by 
the Larmor formula Eq. (1) for orbits of mesoscopic ra- 
dius 2\l- 
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with the London penetration depth A^ given by the usual 
formQ 
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and n s = n the superfluid density. This suggests 
that (i) the transition to superconductivity involves an 
expansion of electronic orbits from radius kp 1 to ra- 
dius 2Al, and (ii) that the superconducting ground state 
can be understood as being composed of all electrons 
(n s = n) residing in orbits of radius 2A^. The latter 
((ii)) is supported by the observation that the angular 



momentum carried by the electrons in the Meissner cur- 
rent with velocity v s in a cylinder of radius R » Xl and 
height h can be written in the two equivalent forms: 

LMeissner = n s (2nRXLh) x (m e v s R) 

= n s {irR 2 h) x (m e v s (2X L )) (5) 

The first form is the conventional description of electrons 
flowing in a surface layer of thickness Xl, hence cross- 
sectional area 2ttRXl, each moving in a circle of radius R. 
The second form describes all electrons in the bulk, hence 
cross-sectional area 7ri? 2 , each moving in a mesoscopic 
orbit of radius 2Xl- 

The proposition that the transition to superconduc- 
tivity involves an expansion of electronic orbits from ra- 
dius kp 1 to radius 2A^ provides a dynamical explana- 
tion of the Meissner effect through the Lorentz force act- 
ing on radially outgoing electrons [H, H|, an effect which 
does not exist in the conventional theory of superconduc- 
tivity. We have argued that the conventional theory is 
untenable § because it cannot provide a dynamical ex- 
planation of the Meissner effect @ and because it leads 
to non-conservation of angular momentum^. 

We showed in Ref. [3j that in the absence of an ap- 
plied magnetic field this orbit expansion in the presence 
of the ionic electric field will generate through the spin- 
orbit interaction a spin current near the surface of the 
superconductor, with carrier velocity given by 
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with n the normal direction to the surface of the su- 
perconductor (Spin Meissner effect). The spin current 
originates in the superposition of real space orbits of ra- 
dius 2Xl for electrons moving with speed Eq. (6), each 
electron thus carrying orbital angular momentum S/20. 

Electronic orbit expansion and the associated radially 
outgoing electron flow lead to a non- homogeneous charge 
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distribution in the superconducting state, with extra neg- 
ative charge near the surface, as proposed within the the- 
ory of hole superconductivity [1, M. New electrodynamic 
equations describe this scenario [Tolfllj]. that yield a defi- 
nite relation between the magnitude of the expelled neg- 
ative charge and resulting electric field in the interior of 
superconductors, and the velocity of carriers in the spin 
current Eq. (6) (Ref. [IT| . Sect. 4). These equations lead 
to the surprising prediction that the density of carriers 
of given spin orientation near the surface will change as 
a function of the magnitude of the applied magnetic field 
and resulting charge current (Ref.pjJ, Eq.(26c) and[L2|). 
leading to a non-zero Knight shift in the superconducting 
state, consistent with experiments fl3l] . 

In addition, we have proposed in earlier work that the 
gap function in superconductors has a gap slope of uni- 
versal sign ( Ref. jl4| . Eq. (19) and Fig. 7) (which mani- 
fests itself experimentally in an observed tunneling asym- 
metry of universal sign[15(), which gives rise to a shift in 
the chemical potential in going from the normal to the 
superconducting state of magnitude determined by the 
slope of the gap function [iff, and that superconductivity 
is driven by kinetic energy lowering [rH [l8j and 'undress- 
ing' of carriers [l^, [2(| from the electron-electron as well 



as from the electron- ion interaction [2 lj. 

In this paper we show that all these phenomena can 
be understood as arising from the existence of a hidden 
'hole core' buried deep inside the Fermi sea of metals in 
the superconducting state. 



pointing outward (parallel to n). The electric field Eq. 
(10) is the radial electric field generated at the surface 
of a cylinder of radius r = 2Xl by the ionic charge den- 
sity pi = \e\n s that compensates the superffuid electronic 
charge density en s {n s = superffuid density). The equal- 
ity between Eqs. (9), (10) and the second term in Eq. 
(7) follows from using the usual relation for the London 
penetration depth Eq. (4). In this interpretation, the 
spin-orbit interaction can be understood as arising from 
the unscreened positive charge density \e\n s acting on 
the macroscopic superfluid wave function. Instead, for a 
single electron in the presence of the electric field 
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that arises from charge expulsion [ll|, the Dirac spin-orbit 
interaction gives the second term in Eq. (7) reduced by 
the factor r q q Q , with r q — h/2m e c the 'quantum electron 
radius' [H. 

The dispersion relation Eq. (7) gives rise to two 
'Rashba-like' bands, given by 
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II. HOLE CORE 

If an electric field exists in the interior of 
superconductors it is only natural to expect that the 
Rashba spin-orbit interaction (2^1 will play an important 
role. We assume that in the normal state electronic en- 
ergies are given by the free electron dispersion relation 
£k = h 2 k 2 /2m e and that in the superconducting state 
spin-split bands develop, with dispersion 
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where h is the outward-pointing normal to the surface of 
the superconductor and 
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The energy dispersion Eq. (7) can be understood as aris- 
ing from the Dirac spin-orbit interaction [23] 



corresponding to spin orientations a ■ (k x h) = ±1. The 
bands are shown in Figure 1. There are negative energy 
states in band 1 in the range < k < qo. We pos- 
tulate that as the system goes superconducting the spin- 
degenerate band splits into the two bands Eqs. (12), (13) 
and simultaneously the negative energy states in band 1 
fill up with holes, driven by the Coulomb interaction [24j| . 
In other words, electrons in these negative energy states 
get ejected from the bulk of the superconductor towards 
the surface. The occupation of the bands in the Brillouin 
zone is shown in Figure 2. The gap between energy states 
in both bands for the same wavevector is 
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as obtained in Ref. [|| using different arguments. 

The carrier density for a two-dimensional free electron 
band structure is given by n s = kp/2ir, the Fermi energy 
by e F = fi 2 kp/2m e , the total energy (per unit area) by 
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with p = ftk and electric field 
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and the energy per particle by e = /2 = h 2 k 2 F /Am e . 
For the spin-split bands, assuming no change in the den- 
sity of particles, the new Fermi wavevectors are 
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FIG. 1: Spin-split bands below T c given by Eqs. (12), (13). 
kF and €f are the Fermi wavevector and Fermi energy (chem- 
ical potential) in the normal state with energy dispersion re- 
lation indicated by the dot-dashed line. In the superconduct- 
ing state the Fermi energy shifts down by v and the Fermi 
wavevectors for the two Rashba bands are £;.fi and The 
negative energy states in band 1 in the region k < qo are filled 
with holes as indicated by the + signs. Note that in the re- 
gion &F2 < k < kFi only electrons of helicity a ■ (k X z) = 1 
exist. 
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and the new Fermi energy is 
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Thus, the Fermi energy (or chemical potential) is lowered 
by v in the transition to superconductivity. This change 
in chemical potential as the system goes superconducting 
was derived within the theory of hole superconductivity 
as originating in the slope of the gap function[16| 
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and coincides with the result Eq. (18) derived here from 
the condition of constant density of particles. 

The energy of electrons in the spin-split bands (per 
unit area) is given by 



2tt2to p 1 4 



'- o / 



(20) 



with a i = 1,<72 = —1. Thus, electrons in band 1 lower 
their kinetic energy upon spin splitting and those in band 
2 raise their kinetic energy. The net change in energy per 
particle is 
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FIG. 2: Brillouin zone occupation for the spin-split bands 
of Fig. 1. The h axis points out of the paper. A circle of 
radius go = l/2Ar, in the (k x ,k y ) plane centered at (0,0) 
is singly occupied by electrons of helicity a ■ (k x n) = — 1. 
The arrows near the center indicate the direction of (k x ,k y ) 
(radial arrow) and of the spin of the corresponding electron 
(attached tangential arrow). There is an excess of electrons of 
helicity (?• (fexn) = 1 at the Fermi surface in the antibonding 
regions close to (k x , k y ) — (±7r, ±7r), (their spin orientation is 
indicated by the arrows). The resulting regions with singly- 
occupied electrons are denoted by the hatched areas, with 45° 
hatching (-45° hatching) in the antibonding (bonding) regions 
respectively. The region hatched by the square pattern is 
doubly occupied with electrons, and the white regions near 
(k x , k y ) = (±7r, ±7r) are doubly occupied with holes. 



in agreement with the result obtained in Ref . [H using en- 
tirely different arguments. Eq. (21) is the condensation 
energy per particle in the superconducting state. Thus, 
the condensation energy in this model can be interpreted 
as arising from kinetic energy lowerin g (in band 1), in 
agreement with previous considerations]!]], [l8| . 
The number density in each band is given by 
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The last term in Eq. (22) is the missing density of carriers 
in the negative energy states, which we argue is expelled 
to the surface of the superconductor 
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assume that the current density is given by 



FIG. 3: Spin-split bands in the presence of a charge cur- 
rent. The occupation in band 1 shifts (the occupied states 
are depicted by the solid line) while that in band 2 remain 
unchanged. The situation depicted corresponds to the critical 
current where the wavevector shift is qo. The position of the 
hole core in the Brillouin zone is shown at the top. 



The speed of carriers in the spin-split bands is given 
by v\ — (l/H)de l k /dk, hence 
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Therefore, the speed of the carriers at the large Fermi 



surface is the same for both bands (v\ Fl = v^, F2 ) 
ever, at the small Fermi surface (k — qo) we have 
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Thus, we can interpret Eq. (27) as the spin current speed 
of the condensate as a whole. In a cylindrical geometry, 
the carriers near the surface with spin orientation paral- 
lel or antiparallel to the cylinder axis will have a 'drift' 
velocity given by parallel to the surface and per- 
pendicular to the cylinder axis, as proposed in our paper 
on the Spin Meissner effect based on entirely different 
arguments [3j. 

When a charge current circulates, we argue that it 
is only the electrons in band 1 that shift their occu- 
pation. This is based on the fact that experiments on 
rotating superconductors [25l |2o| can be interpreted as 
indicating 27] that the superconducting current is car- 
ried by completely 'undressed' electrons, suggesting that 
a shift in the small Fermi surface describing long wave- 
length free-electron-like carriers is involved. Hence we 
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involving only half the superfluid density (n s /2) corre- 
sponding to the carriers in band 1, moving with super- 
fluid velocity u*. The second equality in Eq. (28) follows 
from London's equation. Thus Eq. (28) yields for the 
drift speed of the carriers in band 1 
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corresponding to a shift in their wavevector 
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When the wavevector shift reaches the value qo the small 
Fermi surface crosses the origin and it is reasonable to 
expect that this corresponds to the critical value of mag- 
netic held where superconductivity will be destroyed. Eq. 
(30) then yields for the critical field (when Aq — qo) 
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which is essentially H c i , the lower critical field of type II 
superconductors [2]. The result Eq. (31) coincides with 
the result obtained analyzing the Spin Meissner effect 
based on entirely different arguments [3]- Figure 3 shows 
the occupation of the bands in the presence of a charge 
current. 

In connection with the Spin Meissner effect we 
showed[|J that the spin current with speed Eq. (6) devel- 
ops when carriers expand their orbits from a microscopic 
length to radius 2Al and are deflected azimuthally by 
the spin-orbit interaction. This corresponds in the cur- 
rent scenario to the dropping of holes in band 1 from 
the top to the bottom to form the hole core. The two 
representations of the same physics are shown in Fig. 4. 
Furthermore we showed in refs. [3j, |4j that the orbit ex- 
pansion in the presence of a magnetic field gives rise to 
the charge current needed for the Meissner effect, and 
argued that the orbit expansion is driven by kinetic en- 
ergy lowering. Clearly, in the current scenario this cor- 
responds to the carriers in band 1 developing a charge 
current as the holes condense from the top of the band 
to the hole core in the presence of a magnetic field, with 
the kinetic energy of carriers in band 1 being lowered as 
described by Eq. (21). 

The scenario depicted in Fig. 3 also shows that when a 
charge current circulates the magnetization of the system 
will change. Assuming that the ejected electrons from 
the states indicated by ' + ' in Fig. 3 represent the ex- 
cess electrons giving rise to the negative charge density 
p_ predicted to exist near the surface by the modified 
electrodynamic equations Fig. 3 shows that as the 
critical magnetic field is approached the excess electrons 
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FIG. 4: In real space, orbits expand from a microscopic length 
to radius 2\l (upper left). The upper right picture shows the 
electronic orbits in a cylinder viewed along the axis, electrons 
with spin out of (into) the paper circulate clockwise (coun- 
terclockwise) with speed (Eq. (1)) . The lower left panel 
shows pairs of holes of helicity a ■ (k x n) = 1 dropping from 
the top of band 1 to the bottom to form the hole core, shown 
in the lower right panel, with directions for k and a for holes 
in the hole core indicated schematically. 



moving in direction opposite to the charge current all 
have the same spin orientation, and the excess electron 
density of opposite spin orientation has been depleted to 
zero (since all the holes have moved to the branch on the 
right in Fig. 3). This is precisely the situation predicted 
by the electrodynamic analysis of ref.[ll| using entirely 
different arguments. We found in that work that the 
expelled charge density for spin a is 
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so that the expelled charge density for a parallel to B 
goes to zero at the critical field Eq. (31), in agreement 
with Fig. 3. The variation of spin density with mag- 
netic field indicates that the superconducting state has a 
triplet component and willgive rise to a non-zero Knight 
shift at zero temperature [lj]. It is remarkable that the 
electrodynamic analysis of Ref.[ll| predicted this 'Edel- 
stein effect '[28], which we can now understand as origi- 
nating in the locking of spin and momentum variables in 
the Rashba-like bands. 

Finally we would like to show that the expelled elec- 
trons from the negative energy states of band 1 corre- 
spond to the excess negative charge near the surface de- 
rived from the electrodynamic equations, given bv(ll| 
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This excess charge is distributed over a surface layer of 
thickness Al, hence the 'surface charge density' E is 
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The charge content of the hole core (Eq. 24) is 
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with a 



2 /he = 1/137 the fine structure constant. 



Thus, Eq. (35) implies that the expelled negative charge 
density p~ originates in the hole core of 137 lattice planes, 
which itself is of the order of A^, as first suggested by 
Slater [H. 

Alternatively we can argue: if Az is the spacing be- 
tween lattice planes parallel to the surface, the three- 
dimensional charge density resulting from Eq. (24) is 
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and setting n^~ p = p-je yields Az — a\h- Taking for 
the spacing between lattice planes Az = 2ao, with a$ = 
h 2 /m e e 2 the Bohr radius yields 
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corresponding, for a simple cubic lattice of lattice con- 
stant a = 2ao, to a band filling 1/27T = 0.16 holes per 
atom, an almost filled band, consistent with the assump- 
tions of the theory of hole superconductivity. 



III. DISCUSSION 

Within the theory of hole superconductivity, negative 
charge is expelled from the interior to the surface of a 
metal in the transition to superconductivity. A spin cur- 
rent exists in the superconductor near the surface, and 
the superconducting state can be understood as a super- 
position of electronic orbits of radius 2A^. Superconduc- 
tivity is driven by kinetic energy lowering, the chemical 
potential drops in the transition to superconductivity by 
an amount v related to the slope of the gap function ver- 
sus energy which has universal sign, and v also gives the 
condensation energy per carrier. The spin density near 
the surface varies with the applied magnetic field and as- 
sociated superfluid velocity. Carriers 'undress' from the 
electron-ion interaction, expand their wavelength and be- 
come free-electron-like as the system becomes supercon- 
ducting. 

All these properties were deduced in earlier work based 
on different arguments. In this paper we have shown that 
they can all be understood within a scenario where a 'hole 
core' of radius qo = 1/2 with hole carriers of helicity 
a-(kxn) — 1 develops at the bottom of the electronic en- 
ergy band. The direction h is the direction of the interior 



6 



+++ + + + + + + + + + + ^+ 
++ + + + + +1+ + + + + + + 

E 

FIG. 5: The plane shown is the boundary of the supercon- 
ductor. An electric field in the interior of the superconductor 
points in the n direction. There is excess negative charge 
within distance Al from the surface of the superconductor 
that was expelled from its interior that screens the interior 
electric field, and there is a spin current near the surface. The 
long and short perpendicular arrows denote the direction of 
the momentum and the dominant carrier's spin. 

electric field pointing towards the closest surface. The 
electric field arises due to negative charge being expelled 
from the hole core states to the surface of the supercon- 
ductor. Thus, the scenario proposed here is consistent 



with and supports the Spin Meissner effect scenario [3j 
and its relation with negative char ge e xpulsion deduced 
from the electrodynamic equations The 'hole core' 
is the k-space representation of this physics. 

Note that in the scenario proposed here, 'Cooper pairs' 
formed by time-reversed carriers, (k — k 4-), are essen- 
tial, with the spin directions t and 4- corresponding to 
a ■ (k x n) — ±1 rather than being fixed in space. Note 
also that the number of holes in the hole core is only 
~ 10 _6 n s , so that the energy cost in promoting those 
electrons to the Fermi energy is only of order fj,eV per 
supcrfluid electron. 

We believe that a full description of the superconduct- 
ing state of matter has to include the physical elements 
discussed here (which are not part of conventional BCS 
theory) and that this physics should be experimentally 
verifiable. In a planar geometry, the direction of the elec- 
tric field in the interior and the direction of the major- 
ity electron spin for given k direction is shown in Fig. 
5. The k-space scenario shown in Figure 2 should be 
testable by hi gh resolution spin-resolved angle-resolved 
photoemission [30l |31| . keeping in mind that depending 
on the particular band structure the location of the spin- 
split regions could be different from those shown in Fig- 
ure 2. 
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